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1. Introduction 

The short distance behaviour of the partons inside strongly polarized nucleons constitutes one of the 
central research topics in QCD being explored both with perturbative and non-perturbative meth- 
ods. During the last years the polarized deep-inelastic scattering data have further improved [1-15]. 
In this note we report on a new next-to-leading order (NLO) QCD analysis of these data [16], 
updating earlier investigations [17]. At large enough four-momentum transfer Q 2 = —q, the 
structure function gi(x,Q 2 ) mainly receives twist-2 contributions 1 and is related to the polarized 
twist-2 partem distribution functions (PDFs). We analyze the structure function g\(x,Q 2 ), which 
is derived from the longitudinal polarization asymmetry accounting for a data-based description 
of the denominator function [19] and corresponding parameterizations for the longitudinal struc- 
ture function, cf. [16]. In the present analysis we include the 0(a s ) contribution due to charm 
quarks, [20,21] 2 . The structure function g2(x, Q 2 ) is described at leading twist by the Wandzura- 
Wilczek relation [21,23]. The parameters of the polarized partem densities, which can be measured 

Nf=A 

using the above data sets, are determined with correlated errors along with the QCD scale Aq^ d . 
We also analyze potential contributions of higher twist. Finally, a phenomenological parameteriza- 
tion of the polarized NLO parton distribution functions is provided in terms of grids for the central 
values and correlated errors, [24]. 

2. The Analysis 

The NLO QCD analysis of the structure function g\ (x, Q 2 ) is performed in Mellin space following 
the standard formalism, cf. e.g. [25], including the heavy quark corrections [26]. 3 In this represen- 
tation the evolution equations can be solved analytically, in both a fast and numerically precise way. 
Only one numerical integral around the singularities of the solution in the complex plane, located 
on the real axis left to an upper bound, has to be performed to represent g\(x, Q 2 ). As the data 
are located at low values of Q 2 target mass corrections are applied, cf. [18, 28]. For the deuteron 
targets a wave function correction is performed [29]. The parton distributions at the starting scale 
Qq = 4 GeV 2 are parameterized by 

xAfi(x,Q 2 ) = ifcV(l -x) bi {\ + Yix) , (2.1) 

with r\i the first moments. The present analysis parameterizes the sea quarks assuming approximate 
flavor St/ (3) symmetry. The deep-inelastic data alone cannot resolve the flavor dependence of the 
sea. Taking into account semi-inclusive data [30], and later on polarized Drell-Yan and di-muon 
data, will allow the determination of polarized sea quark distributions, similar to the unpolarized 
case. 4 r\ Uv and 7}^ are fixed due to the neutron and hyperon-/? decay parameters F and D, which 
are very well measured : 

H«, ~rid v =F + D and r) Uv + r) dv = 3F + D, (2.2) 



'Twist-3 contributions are connected by target mass effects, cf. [18]. 

The O(a^) heavy flavor corrections are only known in the asymptotic region Q 3> m , [22]. 
3 We refrain from carrying out small-x resummations, since yet unknown subleading terms are very likely to cancel 
the leading order effects, cf. [27]. 

For a first analysis accounting for the flavor dependence of the sea quarks see [31]. 
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T7„ r =0.928 ±0.014 and r\ dv = -0.342 ±0.018 . 



(2.3) 



The parameters in (24) cannot all be measured using the present data since for some the % 2 - 
fit yields errors larger than 100%. In case of the sea-quark and gluon density ji is found to be 
compatible with zero. Furthermore the ^-parameters of the distributions Aq s and AG are related by 
about a ag = aA q , + 1, which we use. The parameters y Uv and y Uv are fitted in an intial run and are 
then kept fixed as model parameters. For the large-* parameters b^ and b^c we used the relation 
£>a ?s /^ag(P 01 ) = i>Aq s /^AG( un P°l) = 1-44 and determine b^c = 5.61 and b^q, = 8.08 in the fit. In the 
final fit 8 parameters are determined including Aq^ . In Figure 1 we show the four distributions 
Au v ,Ad v , Aq s and AG at the input scale and compare them to other determinations. In Ref. [16] we 
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Figure 1: NLO polarized parton distributions at the input scale Qq = 4.0 GeV 2 (solid line) compared 
to results obtained by GRSV (dashed-dotted line) [32], AAC (dashed line) [33], and LSS (long dashed 
line) [34]. The shaded areas represent the fully correlated 1(7 error bands calculated by Gaussian error 
propagation. The dotted line indicates the positivity bound using the parameterization [35]; from Ref. [16]. 

provide the correlated errors. Due to this one may perform Gaussian error propagation for all 
observales based on polarized parton denstities predicting the PDF errors of these quantites. There 
we also compute a series of moments for the different parton densities which can be compared to 
upcoming lattice simulations. 

The nucleon spin is given by the relation 



^ = ~{AL(x)) + (AG(x))o + L q + L g 



(2-4) 
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Figure 2: A summary of the current measurements of a s (M^) from unpolarized and polarized DIS data, 
cf. Ref. [16] for details. Due to the size of errors we include only the results of NNLO and N 3 LO analyses 
in the unpolarized case, while those in the polarized case stem from NLO analyses. The yellow band marks 
the world average of 2009 [36]. 

to the first moments of the polarized flavor singlet and gluon distributions and the quark and gluon 
angular momenta L qg . In the present analysis we obtain 



(AZ(x)) 
(AG(x)) 



0.216 ±0.079 
0.462 ±0.430 



(2.5) 
(2.6) 



which saturates the required value even for vanishing values for L q and L g . However, the error on 
the gluon density is still rather large. 

In deep-inelastic QCD analyses it is important to determine the PDF-parameters at the initial 
scale <2q together with the QCD scale Aqcd since there are strong correlations, e.g. between the 
gluon-normalization and a v (M|), but also to other parameters. We obtain 



A 



(4) 
QCD 



243 ±62 (exp) + ^ (FS) ^ (RS) MeV . 



(2.7) 



The renormalization (RS) and factorization scales (FS) were varied by a factor of 2 around Q 2 . 
Here we excluded values /x? < 1 GeV 2 , unlike in Ref. [17], since at scales lower than 1 GeV 2 the 
perturbative description cannot be considered reliable anymore. Correspondingly, one obtains 



«.v(M|) 



^,,™ ±0.0043 . . -0.0029 _ , ±0.0032 _„, 

0.1132 (exp) (FS) (RS) 

-0.0051 v F; ±0.0015 v ; -0.0075 v ; 



(2.8) 



The eiTors are much larger than in the unpolarized case, at NNLO, where an accuracy of 0(1%) 
is reached. Still the central value is lower than the current world average and well comparable to 
the unpolarized values. In Figure 2 we summarize the current status of a s (M^) measurements in 
deep-inelastic scattering. We would like to mention the results of the unpolarized NS-analysis [37] 
at N 3 LO : 

±0.0020 
-0.0022 



CC S (M. 



0.1141 



(2.9) 



QCD Analysis of the Polarized World Data Johannes Bliimlein 



and recent combined NS and singlet NNLO analyses [38, 39] 

a,(Mf) = 0.1124±0.0020 (2.10) 

a s (MJ) = 0.1135 ±0.0014 . (2.11) 

Very recently, due to the inclusion of the combined Hl+ZEUS data the latter value receives a slight 
change to 

a s (Ml) =0.1147 ±0.0012 , (2.12) 

reaching now the accuracy of 1 %, [40]. 

We also fitted additive higher twist terms to g\ (x, Q 2 ) to explore the corresponding structures 
in the region x < 0.6 for the proton- and deuteron targets. While in case of the deuteron target the 
results is fully compatible with zero, an effect of up to 2 a is observed for four of five bins in case 
of the proton target. This result is indicative mainly. In the unpolarized case one has much better 
means to separate leading and higher twist effects, cf. [41,42], in a clear manner. This also requires 
even higher order corrections. A comparable analysis in the polarized case has to be based on much 
more precise data in a far wider range of Q 2 which can be obtained at future colliders such as the 
EIC. 
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